Effects of measurement back-action in the stabilization of a Bose-Einstein condensate 

through feedback 
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We apply quantum filtering and control to a particle in a harmonic trap under continuous position 
measurement, and show that a simple static feedback law can be used to cool the system. The final 
steady state is Gaussian and dependent on the feedback strength and coupling between the system 
and probe. In the limit of weak coupling this final state becomes the ground state. An earlier model 
by Haine et. al. (PRA, 69 2004) without measurement back-action showed dark states: states that 
did not display error signals, thus remaining unaffected by the control. This paper shows that for a 
realistic measurement process this is not true, which indicates that a Bose-Einstein condensate may 
be driven towards the ground state from any arbitrary initial state. 

PACS numbers: 03.75.Hh,02.30.Yy,02.50.-r,42.50.Lc 



I. INTRODUCTION 

The observation and control of a Bosc-Einstcin conden- 
sate (BEC) using a real-time feedback loop has become 
an important scientific objective within the atom optics 
community. Such a process may be essential to reduce 
the noise characteristics of an atom laser that is outcou- 
pled from such a BEC. There are two classes of noise that 
require control. The first is the quantum statistical noise 
of the BEC mode, which limits the lincwidth of an atom 
laser outcoupled from a sing lc mode BEC [1 0. The 
second is the spatial mode of the BEC itself, which may 
not be stable under varying preparation conditions Q , or 
continuous pumping 0, [j[ • The quantum statistical noise 
of the mode will not be the dominant effect on the prop- 
erties of the atom laser unless the BEC can be produced 
in a single spatial mode, so the first consideration of im- 
plementing control schemes will be to provide that stable 
spatial mode. This paper shows that a straightforward 
static measurement feedback can address that problem. 

In previous work, Haine et al. address controlling a 
BEC system with and without nonlinearities using re- 
alistic trap controls, but did not model the back-action 
of their measurement system They model the sys- 
tem semiclassically, and provide feedback based on the 
state of the system, rather than include the effects of cou- 
pling the BEC to a measurement device. In the absence 
of atomic interactions, this work predicted non-Gaussian 
'dark states' that did not respond to the control, although 
they were not in the ground state. Doherty and Jacobs 
include this back-action in a model of the measurement 
of a single atom trapped in a harmonic potential, and 
they provide an optimal control scheme for driving it to- 
wards the ground state [f|. This result only holds for a 
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Gaussian initial state, so even in the absence of atomic 
interactions, the dark states of the semiclassical treat- 
ment could not be addressed by that analysis. Position 
measurement and control schemes for Gaussian quantum 
systems have also been considered for nanomechanical 
resonators @ and atomic cavity QED 0]. 

In this paper we show that a single trapped parti- 
cle in an arbitrary quantum state can always be cooled 
down to the ground state by the experimentally realistic 
feedback scheme shown in reference Q, as measurement 
back-action destabilises the dark states shown in that 
analysis. 

In Sec. [Ill we introduce our model of the system and 
the quantum filter that is used to provide the control. 
Sec. IIIII produces analytic results for the final expected 
energy of the system under the influence of feedback. It 
does this first for a Gaussian initial state, and then shows 
that the final state is independent of the statistics of the 
initial state. Sec. IIVI shows the details of the numerical 
simulations that demonstrate this process. 



II. MODEL 

Modelling the full quantum field of a BEC in the pres- 
ence of interactions between the particles leads to a model 
that is extremely difficult to solve either analytically or 
numerically. Fortunately, the pathological case in the 
semiclassical model was not the interacting case, but the 
linear case [4j| . In the absence of interactions there were 
non-stationary states that were not cooled by the feed- 
back. By contrast, the presence of interactions coupled 
these states to others that were affected by the feedback, 
and the system approached the ground state. We inves- 
tigate the linear case therefore, where the BEC is equiv- 
alent to a single trapped particle, and the full quantum 
system is tractable. 

The trapping potential for the BEC is modelled as a 
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harmonic potential of angular frequency w, with a con- 
trollable linear potential u(t)q. In the absence of inter- 
actions, the Hamiltonian is given by: 
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where u(i) is the control signal, to be determined. As 
the Hamiltonian is quadratic in system operators and 
bilinear in system operators and the control signal, it is 
expressed in the more compact matrix form 
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where G is real and symmetric, B is real and x is a vector 
of system operators such that the commutation relations 
are represented by xx T — (xx T ) T = ihH. For the system 
considered here these become 
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The feedback control of this system via a process 
of continuous measurement requires a model of the in- 
teraction with the measurement device. There are 
many treatments in the literature of how to introduce 
a measurement on a quantum s yste m in different con- 
texts H, ©, 01, El El El Q 111- We consider the 
interaction between the system and the probe field to 
be represented by the coupling operator L = L T x, such 
that the master equation for the continuously monitored 
system will then be 

dp = ^[H A ,p] dt + V[L]pdt 
n 

= ^[H A ,p] dt + aV[q]pdt, (4) 
n 

where T>[c]p = cpc< — \j2(c) cp + pc^c) and the definition 







for L, where a is a real-valued number, has been used 
in the second line. This represents a system undergo- 
ing a continuous position measurement, which can be 
obtained, for example, by placing a BEC in an optical 
cavity ||. Here a is the measurement strength and de- 
pendent on physical parameters of the probing system 
(e.g., optical cavity) but is left unspecified. It is also 
formally possible to extract position measurements from 
system-reservoir couplings that equate to pure damping 
of the motional state, but it is precisely the absence of 
such cooling mechanisms that motivates our search for 
control schemes. This position measurement produces 
momentum diffusion in the system, and we will show 
that this momentum diffusion makes the ground state 
an attractor for the system under feedback control. 



Measurement of the real quadrature of the outgoing 
probe field provides a measurement signal dq meas upon 
which the conditional state 7r t is updated. Thus the con- 
ditional density operator Tr t is a function of the mea- 
surement signal, and provides a best estimate X t of any 
system operator X conditioned on those measurement 
results (in the mean square sense): 

X t = Tr[n t X}. 

The ensemble average E[Aj] is equal to the expectation 
calculated from the entire model (X t ): 



E[X t ] = Tr[p Q X t ] = (X t 



(5) 



so that E[7r t ] — p t i.e. the classical average of the condi- 
tional density operator 7r t is equal to the unconditional 
density operator p t . Using either the reference probabil- 
ity techni que fill [ItJ , or the innovations method [H, [l9| 
(see also [20, [21] ) , the stochastic master equation for the 
quantum filter is 

— i 

d-K t = — [H A ,Trt] dt + aV[q}Tr t dt 

+^Jr\a (qiTt + TT t q - 27r 4 Tr [qn t ]) dW, (6) 

where dW/2^frja = dq meas - q t dt = dq meas - Tr [ir t q] dt, 
dW is gaussian white noise (a Wiener process) and 77 is 
the detection efficiency. The stochastic term in the above 
equation represents the additional noise introduced by 
the measurement signal dq meas due to the finite strength 
of the continuous measurement. 

A control signal u(t) is fed back into the position of 
the minimum of the trap potential, and this involves a 
term in the Hamiltonian of the form H contr oi = u(t)q. 
Note that for arbitrary system operator X the dynamical 
equation for the conditional expectation X t has a term 
due to the control 
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TTt] X 



u(t)Tr[[X,q]n t ]. (7) 



Thus a control proportional to q affects an operator X 
only if the two do not commute, and controlling the trap 
position via the term u(t) q involves feedback to the mo- 
mentum quadrature. 

We now define the control law u(t) as a simple feedback 
of the conditional expectation of position and moment 
qt,Pt in a linear form: 

u(t) = k q q t + kpp t . 

This is s state-estimate feedback, in line with the choice 
used by [3| who feed back the position moment, and is 
similar to that chosen by [BJ except that the requirement 
of optimality has been relaxed. Using the vector notation 
of Eq.Q it is fruitful to define a vector of conditional 
means 



x t = Tr [x TT t ] , 



(8) 
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so that the control term in the Hamiltonian can be ex- 
pressed as 



Bu 







) 


( «(*) ) 




k p J 
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(9) 



III. FINAL EXPECTED ENERGY 

The objective for this control scheme is to drive the 
energy of the system towards the ground state, and the 
effectiveness of the control can be discerned from the 
final energy of the system. Since this is a quantum 
stochastic formulation, the expectation, or averaging, 
over all possible quantum trajectories corresponding to 
all possible measurement outcomes is required. More- 
over for a general state (other than an energy eigen- 
state) the final energy will be a random variable that 
will require an expectation. This quantity can calculated 
from the energy operator in the Schrodinger picture given 
by E — ix T Gx and the quantum filter 7r t , using the 

property (E t ) = E[E t ] = E [Tr [n t , E]]. The long term 
stochastic limit for the expected energy is of interest, 
which requires the stochastic master equation for the fil- 
ter Eq.Q and the equation for the feedback term given 
by Eq.©. This limit for a Gaussian state is analytically 
calculated below, and an argument that all initial states 
will be driven towards this same limit is provided, along 
with supporting numerical simulations. 

Insight is gained by looking at the evolution of a Gaus- 
sian density operator under Eq.([S]). If an initial state is 
Gaussian then it remains Gaussian, due to the linearity 
of the coupling, the quadratic form of the Hamiltonian 
and the quantum Gaussian noise. 

Assuming a Gaussian form for the density operator 
TTt , it can be described by a vector of means x t given by 
Eq. ([5]) and a matrix of symmetrized covariances V t given 
as 



V t = iTr[(xx 2 
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X t X t . 
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Using Eq.© and the control law Eq.@ the dynamical 
equations for the means and symmetrized covariances are 
calculated in their matrix form and come out as 

dx t = {A + K)x t dt + ^ (2V t L) dW (11) 
V t = AV t +V 4 A T + D-4ryV 4 L(L) T V t , (12) 

where D = ft 2 S [L(L) T ] S T and A = EG, and 
EqfT2l is often referred to as a Riccati equation. Using 
the definition of the covariance matrix ~V t Eq. lfTQ]) the 
average of the conditional energy provides the expected 
energy of 

(Ek) = ^E [xf Gx t ] + iTr [GV t ] . (13) 

The expected energy is proportional to the expectation 
of the square of both position and momentum operators, 



which, for a Gaussian operator, is given by the variance 
plus the means squared. 

It is possible to analytically solve Eq. (fT3")) in the limit 
t — ► oo. Eq. ljlip for the means is a multivariate Ornstein 
Uhlenbeck process with time dependent parameters [13] • 
Note that the variance matrix V t changes with time. For 
a negative value of k p the eigenvalues of (A + K) has 
strictly negative real parts, meaning it is a stable or Hur- 
witz matrix. This implies the dynamical system Eq. (jll|) 
is stable and exp[(A + K)i] — > as t — > oo. The solution 
is formally written as 

x t = exp [(A + K)t] x 

+2^j f exp[(A + K)(i-s)]V a LdW..( 14 ) 
Jo 

Since having k q ^ does not provide any damping of 
the means, for the rest of the paper it is assumed k q = 
and k p < 0. This is easily interpreted: a control signal of 
a negative momentum estimate, coupled to the position 
operator, feeds back into the momentum estimate, and 
thus the control provides damping of the momentum, and 
both of the conditional means. 

Convergence of the Riccati equation Eq. (fT2|) to a 
unique non-negative definite stabilizing solution is as- 
sured in our case [23j, and can be obtained by solving 
the algebraic Riccati equation, Voo = 0: 
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It is fortunate that this equation has an analytical solu- 
tion 0, 
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where the steady state for variances are given as @ 
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It is easy to check these variances obey the Heisenberg 
Uncertainty principle, and the purity of the state can be 
calculated as 0, Hj] 
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= Vv- 

Expanding out the variances to 2nd order for small a 
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Thus the steady state limit t — > oo of the means and 
variances are available. Since the Gaussian state can 
be described by these parameters and their limits, the 
stochastic steady state limit under Eq. ([6]) for a Gaussian 
state is possible. However, using just the long term limits 
of both the variances Eq. (fT2|) and the expectation of the 
square of the conditional means Eq. (fT4")l [22] the limit for 
expected energy Eq.(TT3|) comes out analytically as 

(£oo> = \mu J 2 V^ + VS>/2m 
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As (-E00) has terms that are both linear and inversely 
proportional to the feedback strength k p , there exists an 
optimal negative value for the strength of the feedback, 
as shown in Figfl] and given by 
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For all the numerical simulations shown in this paper 
we used h^m = LO = rj=l and therefore all quantities 
displayed in the figures have dimensions compatible with 
this choice. Fig[2]shows the relationship between the final 
expected energy and the coupling strength between the 
probe and system. It highlights the fact that a weaker 
coupling introduces a smaller disturbance to the system, 
which can also be seen for the expansion to 2nd order for 
small a of (-Eoo) 



(Eoo)(a) 
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If the initial state ttq is not Gaussian, then Eq. (fT5|) 
can still be expected to hold. While a detailed proof of 
this is beyond the scope of the present paper, evidence 
supporting this is given in the simulations in the following 
section. In fact, the limit conditional state is expected to 
be Gaussian, with zero mean and covariance V^. Such 
limits were investigated in [251 ] . and in the classical case 
in 0. 

To see why this might be true, we expand an arbitrary 
initial state ttq as a sum of Gaussian states, i.e. a P 
representation over the coherent states |/3): 



ttq = / P(/3) \f3)((3\df3. 



(17) 
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FIG. 1: Final expected energy, Eq. (|15[l . as a function of the 
feedback parameter k p , for h = m = uj — r] — 1 and a = 0.09. 
The horizontal line represents the ground state energy, while 
the cross represents the calculated optimal control strength. 
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FIG. 2: Final expected energy, Eq. (|15)) . as a function of cou- 
pling strength a, for k p = —1.35. The horizontal line repre- 
sents the ground state energy. 



The stochastic master equation for a state Eq. ([6]) is not 
linear in ir t , so it is not possible to apply the stochastic 
master equation to each state individually. Fortunately 
there exists an equation for the evolution of the unnor- 
malized state at that is linear in a. This is the equation 
for the unnormalized Belavkin filter (l6| . the quantum 
counterpart for the Zakai equation from classical filter- 
ing, 



da t 



—i 11 

~j~ i H A, <Tt] dt + a(qa t q - -<? 2 0( - ^tQ 2 ) dt 

(18) 



-2^/rja (qa t + a t q) dq„ 



where dq meas is a measurement signal driving the system, 
and (To = ?u> In quantum filtering this equation is often 
derived first, and the stochastic master equation Eq.© 
for the normalised state lit is easily achieved through a 
process of normalization: here -Kt — <Jt/nt, where rit = 
TrH. 

Let <rf be evolution of the individual coherent initial 
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state ((3\ = tTg under Eq[l8l with normalisation factor 
rif = Tr erf and corresponding normalised state 7rf = 

crf/nf. By linearity the unnormalised states erf can be 
summed, and the total then normalised, to end up with 
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n, 



(19) 



IV. SIMULATION 



When simulating the system, the stochastic master 
equation for the quantum filter Eq.© is expressed as 
a stochastic differential equation for a Q function [27| 
with the relation 



where the total normalisation factor is given by n t 



Tr 



SP{p)o*dp\ = J P{P)n p t dfi. 
Using the classical limiting results for gaussian initial 
conditions [26], for large times 7rf is independent of (3 i.e. 
?rf sw 7rf 00 for all f3, where 7rf °° is the limiting state given 
by a gaussian of mean x t and variances ~Voo- Thus for 
large times we have 



dQ(x,y,t) = - (x + iy\dn t \x + iy), 

7T 
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(20) where \x + iy) is a coherent state, so that 
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Q(x,y,t) dt 



(21) 



The first term, proportional to ui, causes a rotation 
in phase space around the origin due to the harmonic 
potential. The second term is the effect of the control, 
proportional to k p and the mean momentum estimate, 
and drives the Q function towards the origin. The third 
is the backaction of position measurement, which involves 
a broadening of the Gaussian in the momentum quadra- 
ture. The last part is called the innovation term, and 
involves the Wiener increment dW arising from the ran- 
dom update of the best estimate coming from the mea- 
surement data. It causes the function to be stochastically 
shifted around the phase space as well as a narrowing of 
the state in the position quadrature. 

Equation (|21[) was simulated in a Stratonovich form us- 
ing the numerical package XMDS [28] . The initial condi- 
tion was chosen to be a linear combination of two number 
states, 7T = |^o)(M with |^ ) = (|1) + |4»/V2, which 
satisfies the dark state conditions, where in the absence 
of measurement back-action the feedback would have had 
no effect on the energy of the system Q. The combined 
effect of all the terms in Eq. f2Tj) in the dynamics can 
be seen in Fig. [3] where the Q function at various times 
is plotted for k p = —1.35 and a = 0.3. As argued in 
the previous section, as a result of the measurement pro- 
cess itself the initial dark state (Fig. [3]-a) continuously 
turns into a Gaussian, becoming subject to the effects 
of the control, which will cool the system down. This 
cooling dynamics is shown in Fig.[?]for the same param- 



eters as in Fig. [3] For a single realisation of the noise in 
Eq. ([2"Tj) the energy will asymptotically fluctuate around 
the steady state solution of Eq. (fl"5")) (horizontal line). 
The figure also shows the result averaged over 48 reali- 
sations (dashed line) where the energy starts to converge 
to its asymptotic predicted value. A curve with the so- 
lution, Eq. ([T3"|) . for a given Gaussian initial state is also 
shown for comparison (dot-dahsed line): cooling is effec- 
tive since the beginning for a initial Gaussian while it 
takes longer for an initial dark state. 

V. CONCLUSION 

This paper shows that a quantum system under con- 
tinuous measurement, in a harmonic potential and driven 
by simple feedback law, will evolve towards a stable, low- 
energy state, and that this occurs for all initial states. In 
the limit of weak measurement, the steady state is the 
ground state of the system. This can be deduced from the 
idea that all initial states under continuous measurement 
evolve towards a Gaussian state [25[ , [2!| . We provide an 
analytical result for the steady state limit for a Gaussian 
state under continuous measurement and subject to the 
control law suggested here, which we then infer to be 
the steady state for all initial states. This result opens 
a new possibility to control the spatial mode of a BEC 
as it eliminates concern raised by Haine et al. 0] about 
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FIG. 4: Time evolution of the expected energy for the same 
parameters as in Fig. Solid line represents a single tra- 
jectory, dashed line is the result after averaging over 48 re- 
alisations of the noise in Eq. (|2ip . and dot-dashed line is the 
solution for an initial Gaussian state (see Eq. (|13[) ). The hor- 
izontal line corresponds to the steady state solution given by 
Eq. (|15p and the letters to the times at which the Q-functions 
in Fig. [3] are plotted. 
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FIG. 3: Q-functions at different times (corresponding to the 
letters in Fig. Q for the evolution of the dark state iro (a) 
under Eq.fl2TJ, for k p = -1.35 and a = 0.3. 



dark states: energetic states that do not display an error 
signals and hence are not susceptible to control. These 
states are no longer stable under the measurement and 
control scheme suggested here, and thus can be controlled 
and driven towards the ground state. 



The analytic results are supported by numerical sim- 
ulations based on a phase-space representation of the 
quantum state. These simulations are general in that 
they allow any quantum state for the system and a wide 
range of allowable evolution, but are only tractable for 
single mode problems. These numerical simulations can 
provide the basis for extensions of this linear model to 
non-linear systems. 



The stochastic steady state has a minimum energy at 
the optimal feedback strength, and the overall energy of 
the steady state is limited by the measurement strength. 
Setting the measurement strength to a infinitesimal value 
is counter-productive in practice, of course, as the time- 
scale for the convergence to the steady state becomes 
infinitely long. In practice it may be useful to adjust 
the strength of the measurement to provide rapid con- 
vergence and then a low steady-state energy, but this 
will be examined in future work. 
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